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ABSTRACT 

The 609 /iin 3 Pi — > 3 Po line of neutral atomic carbon has been detected in the 
protoplanetary nebulae CRL 618 and CRL 2688. C I appears to be about 7% as abundant 
as CO in CRL 2688, and 70% as abundant as CO in CRL 618. Cl emission arises 
primarily from the slow component of the wind from these objects. The data presented 
here indicate that Cl gradually becomes more abundant as the star leaves the AGB, and 
C I is significantly enhanced relative to CO before the star develops a large H II region. 

Subject headings: circumstellar matter — stars: AGB and post- AGB — stars: evolution 
— stars: individual: (CRL 618, CRL 2688) — stars: mass-loss 
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1. Introduction 

Shortly after the birth of mm wavelength astron- 
omy, observations of CO in planetary nebulae (PNe) 
showed that these objects often have a neutral enve- 
lope as massive or more massive than the optically 
visible nebula (Mufson et al. 1975, Huggins et al. 
1996). The arrival of large submillimcter telescopes 
ten years ago allowed another neutral gas component, 
atomic carbon, to be studied via its 609/i ground- 
state fine structure line Cl(l — 0). Because the fine 
structure lines of CI arise from transitions between 
three rather closely spaced levels, the partition func- 
tion is relatively insensitive to temperature through- 
out a very large range of plausible temperatures for 
neutral interstellar material. This feature makes is 
possible to estimate the abundance of CI using only 
the strength of the C 1(1 — 0) transition. 

Bachiller et al. (1994) observed Cl(l - 0) in the 
Ring Nebula, and found that atomic carbon was sev- 
eral times more abundant than CO in that object. 
Similarly, Young et al. (1997a) found CI to be more 
abundant than CO in the Helix Nebula. On the other 
hand, the Cl(l — 0) transition is barely detectable 
with current receivers in the prototypical mass-losing 
carbon star IRC+10216 (Kccnc et al. 1993). Because 
AGB stars such as IRC+10216 arc thought to be the 
progenators of planetary nebulae, it is natural to ask 
at what point in the evolution of such an object does 
Cl become a major constituent in the neutral enve- 
lope. The transition from the AGB to the PNe stage 
is thought to take only ~ 1000 years (Iben and Ren- 
zini 1983), so transition objects, or protoplanetary 
nebulae (PPNe), are quite rare. CRL 2688 (the Egg 
Nebula) and CRL 618 are two objects which are be- 
lieved to be PPNe. 

In the mm regime, CRL 2688 and CRL 618 superfi- 
cially resemble AGB stars. Both have thick molecular 
envelopes, exhibiting strong, optically thick CO emis- 
sion. However in addition to the slow, ~ 20 kms -1 
wind typical of an AGB star, these objects also have 
a fast ~ 100 kms -1 wind not seen in the envelopes 
of typical AGB stars (Gammie et al. 1989, Young 
et al. 1992). In the visible regime, the photospheric 
emission from these objects can be seen reflected from 
dust in the nebula's dense inner region. Both of these 
objects have much hotter photospheres than would 
be found on the AGB. CRL 2688 appears to be the 
less-evolved object, because its spectral type is F5 la 
(Crampton et al. 1975) while that of CRL 618 is a 



hotter B0 (Westbrook et al. 1975). 

In this letter I present C 1(1-0) spectra of CRL 2688 
and CRL 618, and compare the results with earlier 
investigations of C I and CO in AGB stars and plan- 
etary nebulae. 

2. Observations and Results 

The observations were made using the 10.4 m 
telescope of the Caltech Submillimeter Observatory 
(CSO) on Mauna Kea, Hawaii, during 1993 Novem- 
ber 8-11. A double sideband SIS receiver was used in 
combination with an AOS spectrometer, which had 
a total bandwidth of 560 MHz and a resolution of 
1.6 MHz, corresponding to 1 kms -1 at the frequency 
of the 3 Pi -> 3 P transition (492.1607 GHz). The 
observations of CRL 618 were made in moderately 
good weather, during which the single sideband sys- 
tem temperature (Tgge) was ~ 5800 K, on aver- 
age. CRL 2688 was observed during periods of ex- 
tremely good weather, during which the (Tssb) aver- 
aged ~ 1900 K. After every 20 seconds of on-source 
integration, one of two reference positions, located 
±180" from the source in azimuth, was observed for 
an equal time. In order to minimize any fixed pat- 
tern noise in the spectrometer, a pattern of frequency 
offsets separated by 5 MHz intervals was applied to 
the nominal frequency. Dithering the frequency in 
this way had the additional benefit of removing the 
sideband ambiguity for any line in the composite spec- 
trum. The data were calibrated using the standard 
hot-sky chopper method, and corrected for the small 
(3%) error introduced by the difference in sky opacity 
of the two receiver sidebands. Saturn was observed to 
determine the telescope's main beam efficiency, 46%. 
All temperatures reported here have been corrected 
for this efficiency. 

CI was clearly detected in both objects (see fig- 
ure 1). In addition, there appears to be another spec- 
tral line visible in the CRL 618 spectrum, and perhaps 
very weakly in the CRL 2688 spectrum as well. The 
line's frequency, measured from the CRL 618 profile, 
is 491.973 ± 0.005 GHz. The microwave spectroscopy 
catalogs of Pickett et al. (1996) and Lovas (1992) were 
searched for possible identifications for this line. The 
only plausible line within ±0.005 GHz in either cata- 
log is H 2 CO(7i 7 -6i 6 ) (491.968936 GHz), 106 K above 
the ground state. Lower transisions of H2CO have 
previously been detected in CRL 618 (Cernicharo et 
al. 1989). Lines ~ 100 K above the ground state can 
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Table 1 

609/x Observations of CRL 618 and CRL 2688 



Source 


a 


5 


Vlsr 


Time a 




r LVW 


T? d 
r HVW 


r H 2 CO 




(1950) 


(1950) 


(kms- 1 ) 


(min) 


(K) 




(K kms" 1 ) 




CRL 618 


04:39:34.03 


36:01:15.9 


-22.3 


167 


0.10 


17.8 


16 


8.6 


CRL 2688 


21:00:19.85 


36:29:44.0 


-33.3 


179 


0.05 


5.1 




0.8 



a Total on-source integration time 

b RMS noise level in blank areas of the spectrum after removal of a linear baseline 

c Line temperature integrated over the region where low velocity CO emission is seen 

d This is the additional flux seen in the "wings" of the CRL 618 spectrum, which may be associated 
with the high velocity emission seen in CO spectra. Note that the total velocity extent of the high 
velocity CO emission is too large to fit within the spectrometer's 560 MHz total bandwidth. The flux 
reported here is merely the integrated flux from -160 to 160 kms^ 1 , minus the flux from the low velocity 
wind region. This value is extremely uncertain, because it would be strongly affected by any curvature 
in the spectrometer's baseline. 

This is the flux integrated over the velocity range of the low velocity CO emission, centered on the 
position where the H 2 CO(7i7 — 6i6) line would fall in the spectrum. 



be expected to appear in CRL 618, which has an exci- 
tation temperature of about 90 K in the region where 
molecular emission arises (Bujarrabal et al. 1988). 
The detection of the H 2 CO(7i7 ~ 6i6) line should be 
regarded as tentative, however, because it lies near 
the edge of the spectrum, where the AOS can show 
spurious features due to the roll-off of its sensitivity. 
Table 1 lists the source coordinates that were used, 
the on-source integration time, the velocity integrated 
intensity of the Cl(l — 0) and H 2 CO(7i 7 — 6i 6 ) lines, 
and the RMS noise level of the spectra. 

Observations of CO emission lines from CRL 618 
and CRL 2688 reveal that both objects have high ve- 
locity outflows in addition to the slow winds typical 
of AGB stars (Gammie et al. 1989, Cernicharo et al. 
1989, Young et al. 1992). For both these stars, the 
velocity extent of the CI emission is nearly identi- 
cal to that of the slow wind. However the "baseline" 
of the CRL 618 spectrum shows a gentle curvature 
which could be the result of low-level emission from 
the high velocity wind. The extent of the high veloc- 
ity wind seen in CO (3 — 2), 380 kms -1 (Gammie et 
al. 1989), is too large to be entirely contained within 



the 560 MHz coverage of the spectrometer used for 
these observations, so these spectra cannot be used 
to accurately calculate the amount of C I in the high 
velocity wind material. 

3. Discussion 

To determine the fraction of carbon in the en- 
velope which is atomic, it is useful to compare the 
C I emission with emission from CO molecules, be- 
cause CO has the largest abundance of any carbon- 
bearing molecule. Such a comparison is hampered, 
however, by the fact that low J rotational transitions 
of CO have a large optical depth in both CRL 618 
and CRL 2688. One is therefore forced to use spec- 
tra from a less abundant species, such as 13 CO, and 
then to divide the result by the fractional abundance 
of 13 C. 

Bujarrabal et al. (1994) present 13 CO(2- 1) spec- 
tra of both CRL 618 and CRL 2688 taken at the 30m 
IRAM telescope. They derive an excitation temper- 
ature (T ox ) for several molecular species, and cite a 
number of references in the literature where estimates 
of T ox have been made. While there is considerable 
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Fig. 1. — The CI spectra of CRL 618 and CRL 2688 are shown. In each spectrum the systemic velocity and extent 
of the Low Velocity Wind, derived from CO observations (Gammie et al. 1989, Young et al. 1992), are shown. The 
location of the line tentatively identified as H 2 CO(7i7 — 6i6) is also shown. 
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scatter amongst the estimates, I will follow these au- 
thors, and adopt T ox = 20 K, which falls comfortably 
near the middle of the T ox estimates. For that tem- 
perature, and assuming LTE, the column density of 
CO is given by 

N(CO) = x 5.6 x 10 14 

x r° T mb [ 13 CO(2 - 1)] dV cm" 2 (1) 

where /(X) is the fractional abundance of species X, 
±V is the velocity extent (in km s -1 ) over which 
emission is seen, and T m b is the main beam brightness 
temperature. Similarly, for CI in LTE at 20 K: 

N(CI) = 1.4 x 10 16 

x / T mb [Cl(l-0)] dVcm- 2 (2) 

These simplistic formulae for column densities as- 
sume that the emitting species is optically thin. Com- 
parisons of the 12 CO and 13 CO spectra of CRL 618 
and CRL 2688 imply that 13 CO transitions are not 
optically thick in these objects (Yamamura et al. 
1994, Yamamura et al. 1995), so the error arising 
from the assumption of low optical depth in equation 
1 is probably small when compared with the uncer- 
tainty in the value of /( 12 C)//( 13 C). The rather low 
(< 1 K) peak antenna temperatures for Cl(l — 0) in 
these objects are most easily explained if the C 1(1—0) 
is also optically thin. However it is possible that the 
C 1(1 — 0) emission arises from an ensemble of clumps 
of high optical depth, with a small combined beam 
filling factor. Using equation 2 would then result in 
an underestimate of N(C I). 

There is no general agreement on the value of 
/( 12 C)//( 13 C) in carbon stars. Estimates range from 
8 (Rank et al. 1974) to 40 (Wannier and Linkc 1978). 
I'll adopt the value of 20 obtained by Wannier et al. 
(1991), which was derived using optically thin transi- 
tions of carbon-bearing species in both CRL 618 and 
CRL 2688. 

The Bujarrabal et al. (1994) 13 CO(2 - 1) spec- 
tra are particularly useful, because their 13" beam is 
quite comparable to the CSO's 15" beam at the fre- 
quency of Cl(l — 0). This similarity of beam sizes, 
combined with the fact that interferometer maps sug- 
gest that both CRL 618 and CRL 2688 are slightly 



resolved, at least in CO, by a 15" beam (Hajian et 
al. 1996, Yamamura et al. 1996), allows the calcula- 
tion of the ratio of C I to CO along the line of sight 
to these stars. Entering the velocity integrated line 
intensities from Bujarrabal et al. and the integrated 
intensities from the Cl(l — 0) spectra presented here, 
into equations 1 and 2, I obtain Cl/CO column den- 
sity ratios of 0.7 and 0.07 for CRL 618 and CRL 2688 
respectively. 

Beichman et al. (1983) have previously reported 
the detection of the C 1(1 — 0) transition in CRL 2688, 
using for their observations an InSb bolometer on 
the KAO, which yielded a 150" beam. They derive 
Cl/CO > 5, a result that is clearly inconsistent with 
the value I obtained. It is possible that CI is geo- 
metrically distributed in some way such that much 
more emission is detectable with a large beam. For 
example, if atomic carbon were distributed in a large 
annulus, with a central hole larger than the CSO's 
15" beam, then the C 1(1 — 0) transition could appear 
much brighter with the KAO's large beam. Because 
the CO emission arises from a region comparable in 
size to the CSO's Cl(l — 0) beam (Yamamura et al. 
1996), this would mean that C 1(1 — 0) emission peaks 
in a region significantly displaced from the CO emis- 
sion. Such a situation would be unprecedented - in all 
stars for which maps of C 1(1 - 0) exist (IRC+10216, 
NGC 6720 and NGC 7027) the Cl(l - 0) emission 
peaks at or very near the position of peak CO emis- 
sion. This is true even when C I arises as a product of 
the photodissociation of CO and other carbon bear- 
ing molecules by interstellar UV radiation, as in the 
case of IRC+10216. Also, the similarity of the shape 
and width of the Cl(l — 0) profile to the CO profiles 
argues against dramatically different CO and C I spa- 
tial distributions. In light of the low signal/noise and 
very limited velocity coverage of the Beichman et al. 
Cl(l — 0) spectrum, and considering their nondetec- 
tion of Cl(l — 0) emission from CRL 618 (which ap- 
pears to me to be the brighter object), I suspect that 
they did not detect CRL 2688. However given the 
strength of the emission reported by Beichman et al, 
and the tremendous improvement of submillimetcr re- 
ceivers since 1983, it would not take much telescope 
time to confirm their detection using a small telescope 
such as 1.2 m Mt. Fuji submillimeter telescope, or 
the University of Texas "focal reducer" (Plume and 
Jaffe, 1995), which converts the CSO into aim tele- 
scope. If the Beichman et al. result were confirmed, it 
would mean the C I distribution in CRL 2688 is quite 
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Table 2 

N(CI)/N(CO) in Carbon-Rich Stars 



Source 


N(CI)/N(C0) 


T c ff 


Source 


N(CI)/N(C0) 


T c ff 


IRC+10216 


0.02 a 


1250 b 


NGC 7027 


0.5° 


170,000 f 


CRL 2688 


0.07 


6500 c 


NGC 6720 


10 s 


110,000 f 


CRL 618 


0.7 


30, 000 d 


NGC 7293 


6 h 


90,000' 



a Kccnc et al. 1993, N(CO) from their 13 CO(2 - 1) profile and equation 1 
b Miller, 1970 c Spectral type F5 la (Crampton et al. 1975) 

d Spectral type B0 (Westbrook et al. 1975) e Young et al. 1997b 

f Malkov et al. 1995 9 Bachillcr et al. 1994 

h Young et al. 1997a 1 Mcndcz et al. 1992 



extraordinary. 

Table 2 shows gives the values of N(C I) /N(CO) for 
a set of objects ranging in type from the cool carbon 
star IRC+10216 to the rather old PNe NGC 7293 
(the Helix Nebula). These objects are all carbon- 
rich, and current theories of stellar evolution suggest 
that they are representative of different stages along 
a single evolutionary path. It appears that while C I 
is a very minor constitute of the envelope of a star 
on the AGB, it becomes comparable in importance 
to CO in the PPNe and early PNe stages, and finally 
is much more abundant than CO in a fully matured 
PNe. 

It is interesting to note from table 2 that 
N(Cl)/N(CO) is higher in CRL 2688 than in 
IRC+10216. The calculation of N(C I)/N(CO) is un- 
certain enough that these estimates taken by them- 
selves might not be convincing evidence that C I truly 
is more abundant relative to CO in CRL 2688 than 
in IRC+10216. However the Cl(l - 0) profile of 
IRC+10216 has a dramatically different shape than 
the CO profile (Keene et al. 1993), and CI appears 
to be confined to a pair of thin shells in that object. 
In contrast, the Cl(l - 0) profile of CRL 2688 more 
closely resembles the profile of a CO transition, sug- 
gesting that the optical depth of C 1(1 — 0) is higher in 
this object, and that C I is not confined to a thin shell. 
Since CRL 2688 has not developed a hot enough pho- 
tosphere to produce a significant H II region (Spergel 



et al. 1983), it appears that the photodestruction of 
molecular carbon reservoirs precedes the ionization of 
hydrogen. 

The similarly of the N(Cl)/N(CO) ratios of CRL 
618 and NGC 7027 is somewhat puzzling. NGC 7027 
has a much hotter photosphere, and a much larger H II 
region. In a very young PNe like NGC 7027, much of 
the carbon liberated by molecular photodissociation 
is quickly ionized and C II becomes extremely promi- 
nent (Liu et al 1996). Perhaps the large abundance of 
C I seen in mature PNe is arises when the slow dissi- 
pation of the AGB wind allows the UV radiation from 
the central star to begin destroying CO and other car- 
bon bearing molecules in clumps throughout the en- 
velope. The similarity of the Cl(l — 0) profile in the 
Helix to that of CO (Young et al. 1997a), which is 
known to be clumpy (Huggins et al. 1992), is evidence 
for this interpretation. 

I thank Valentin Bujarrabal, who made several 
helpful suggestions for improvements to this letter, 
and who pointed out some blunders in an earlier ver- 
sion. Many thanks are also due to staff of the CSO 
for their tireless work in supporting observers at that 
facillity. Research at the CSO is supported by Na- 
tional Science Foundation grant AST93-13929. 
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